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Abstract 

We discuss the isospin symmetry breaking (ISB) of the valence- and sea-quark distributions 
between the proton and the neutron in the framework of the chiral quark model. We assume that 
isospin symmetry breaking is the result of mass differences between isospin multiplets and then 
analyze the effects of isospin symmetry breaking on the Gottfried sum rule and the NuTeV anomaly. 
We show that, although both flavor asymmetry in the nucleon sea and the ISB between the proton 
and the neutron can lead to the violation of the Gottfried sum rule, the main contribution is from 
the flavor asymmetry in the framework of the chiral quark model. We also find that the correction 
to the NuTeV anomaly is in an opposite direction, so the NuTeV anomaly cannot be removed by 
isospin symmetry breaking in the chiral quark model. It is remarkable that our results of ISB for 
both valence- and sea-quark distributions are consistent with the Martin-Roberts-Stirling-Thorne 
parametrization of quark distributions. 
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I. INTRODUCTION 



Isospin symmetry was originally introduced to describe almost identical properties of 
strong interaction of the proton and the neutron by turning off their electromagnetic inter- 
action, i.e., their charge information. This symmetry is commonly expected to be a precise 
symmetry , and its breaking is assumed to be negligible in the phenomenological or 
experimental analysis. This is, in general, true, since electromagnetic interactions are weak 
compared with strong interactions. However, it is possible for isospin symmetry breaking 
(ISB) to have important influence on some experiments, especially its effects on the parton 
distributions. Therefore, it is necessary to analyze it carefully. 

The isospin symmetry between the proton and the neutron originates from the SU(2) 
symmetry between u and d quarks, which are isospin doublets with isospin 1 = 1/2 and 
isospin three-components 1/2 and -1/2, respectively. The isospin symmetry at parton 
level indicates that the u (d, u, <i)-quark distribution in the proton is equal to the d (u, d, u)- 
quark distribution in the neutron. Accordingly, the ISBs of both valance-quark and sea- 
quark distributions are defined, respectively, as 
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ISB at the parton level and its possible consequences for several processes were first 
investigated by one of us It was pointed out that both flavor asymmetry in the nucleon 



sea and isospin symmetry breaking between the proton and the neutron can 
violation of the Gottfried sum rule reported by the New Muon Collaboration 



possibility of distinguishing these two effects was also discussed in detail jg]. 

In 2002, the NuTeV Collaboration [?| extracted sin 2 #\y by measuring the ratios of neutral 
current to charged current v and v cross sections on iron targets. The reported sin 2 6*w = 
0.2277 ± 0.0013 (stat) ± 0.0009 (syst) has approximately 3 standard deviations above the 
world average value sin 2 # w = 0.2227± 0.0004 measured in other electroweak processes. This 
remarkable deviation is called the NuTeV anomaly and was discussed in a number of papers 
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from various aspects, including new physics beyond the standard model 8L the nuclear 



effect [9j], nonisoscalar targets |l0|, and strange-antistrange asymmetry llNl3|. Moreover, 
the possible influence of ISB on this measurement was also studied in a series of papers 



201 ] . However, the correction from ISB to the NuTeV anomaly is still not conclusive. 

The Martin-Roberts-Stirling-Thorne (MRST) group 21] provided some evidence to sup- 
port the ISB effects on parton distributions of both valance and sea quarks and included 
ISB in the parametrization based on experimental data. They obtained the ISB of valance 
quarks as 

k(1 -x) 4 x-°- 5 (x- 0.0909), 
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(2) 



where —0.8 < k < +0.65 with a 90% confidence level, and the best fit value is k = —0.2. 
They also obtained the ISB of sea quarks, as can be deduced from Eqs. (28) and (29) in 
Ref. [21 1, 



Su(x) = ku p (x), Sd(x) = kd p (x), 



(3) 



with the best fit value k = 0.08. 

In this paper, we calculate the ISB of the valance- and sea-quark distributions between 
the proton and the neutron in the chiral quark model and discuss some possible effects 
of ISB. We assume that the ISB between the proton and the neutron is entirely from the 
mass difference between isospin multiplets at both hadron and parton levels. 1 In Sec. Ull we 
compute ISB in the chiral quark model, with the constituent-quark-model results as the bare 
constituent-quark-distribution inputs. Then, we calculate the ISB effect on the violation of 
the Gottfried sum rule. In Sec. IHH we discuss the ISB correction to the measurement of the 
weak angle and point out the significant influence on the NuTeV anomaly. In Sec. IIV1 we 
provide summaries of the paper. 



II. ISOSPIN SYMMETRY BREAKING IN THE CHIRAL QUARK MODEL 



The chiral quark model, established by Weinberg 22J and developed by Manohar and 
Georgi 23j, has an apt description of its important degrees of freedom in terms of quarks, 



1 As mass difference between isospin multiplets, especially that between u and d quarks, is not entirely due 
to charge difference, we refer such effect as Isospin Symmetry Breaking (ISB) instead of Charge Symmetry 
Breaking (CSB) as called in some papers. 



3 



gluons, and Goldstone (GS) bosons at momentum scales relating to hadron structure. This 
model is successful in explaining numerous problems, including the violation of the Gottfried 
e from the aspect of flavor asymmetry in the nucleon sea 



sum ru 
sis 



26- 



24 



25| , the proton spin cri 



28| . and the NuTeV anomaly resulting from the strange-antistrange asymmetry 13], 



and has been widely recognized as an effective theory of QCD at the low-energy scale. 

In the chiral quark model, the minor effects of the internal gluons are negligible. The 
valence quarks contained in the nucleon fluctuate into quarks plus GS bosons, which spon- 
taneously break chiral symmetry. Then, the effective interaction Lagrangian is 



if) {iD^ + V„) + igx^A^l^ + 



(4) 



where 



d 



(5) 



is the quark field and = + igG^ is the gauge- covariant derivative of QCD. stands 
for the gluon field, g stands for the strong coupling constant, and g& stands for the axial- 
vector coupling constant. and are the vector and the axial- vector currents, which are 
defined as 



A„ 



1 



(6) 



where £ = exp(iU/ f), and II has the form: 
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Expanding V M and A^ in powers of II//, one gets — + 0(11/ f) 2 and A^ = idfjl/ f + 
0(11/ f) 2 . The pseudoscalar decay constant is / ~ 93 MeV. Thus, the effective interaction 
Lagrangian between GS bosons and quarks in the leading order becomes 



9a 
' f 



(8) 



Based on the time-ordered perturbative theory in the infinite momentum frame, all particles 
are on-mass-shell, and the factorization of the subprocess is automatic, so we can express the 
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quark distributions inside a nucleon as a convolution of a constituent-quark distribution in 
a nucleon and the structure functions of a constituent quark. Since the rj is relatively heavy, 
we neglect the minor contribution from its suppressed fluctuation in this paper. Then, the 
light-front Fock decompositions of constituent-quark wave functions are 

\U) = yfz~ u \u ) + a 7T+ \d7r + ) + ^|wr°) + a K+ \sK + ), (9) 

v2 

\D) = y/Z~ d \d ) + a„-\uir-) + ^\dn°) + a K o\sK°), (10) 

v 2 

where Z u and Z d are the renormalization constants for the bare constituent u quark \uq) and 
d quark \do), respectively, and |ao,| 2 (a = ir,K) are the probabilities to find GS bosons in 
the dressed constituent-quark states \U) and \D). In the chiral quark model, the fluctuation 
of a__bare constituent quark into a GS boson and a recoil bare constituent quark is given 

as 



29] 



Qj( x ) = [ — p j a /i(y)Qi ( - ) > (11) 

j x y v y j 

where Pj a /i{y) is the splitting function, which gives the probability of finding a constituent 
quark j carrying the light-cone momentum fraction y together with a spectator GS boson a, 

rrii, rrij, and m a are the masses of the i- and j-constituent quarks and the pseudoscalar 
meson a, respectively, and fa = (m^ + rrij)/2 is the average mass of the constituent quarks. 
Mja = ( m i + ^t) /y + ( m a + ^t) / (1 — y) ls ^ ne sc L uare °f the invariant mass of the final 
state. We can also write the internal structure of GS bosons in the following form 

Qk(x) = [ ——V k/a ( — J Paj/i ( — J qi (y 2 ) , (13) 



yi yi \y\j \yi, 

where Vk/ a (x) is the quark k distribution function in a and satisfies the normalization 
Jq V k/a (x)dx = I. 

When we take ISB into consideration, the renormalization constant Z should take the 
form 

Z u = l-(P n+ )-^(P uw0 )-(P K+ ), 

Z d = l- (P w _> - i(P^,) - (P K0 ), (14) 
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where (P a ) = (Pj a /i) = (Paj/i) — Jo xU 1 Pja/i{ x )^ x [29]. It is conventional to specify the 
momentum cutoff function at the quark- GS-boson vertex as 



9 A ->■ 9a gx P 



m 2 j - Mf a i 
4A 2 



(15) 



where g' A = 1, following the large N c argument 30|], and A is the cutoff parameter, which is 
determined by the experimental data of the Gottfried sum and the constituent-quark-mass 
inputs for the pion. Such a form factor has the correct t- and w-channel symmetry, Pj n /Ay) = 
P a j/i{l ~ y)- Then, one can obtain the quark-distribution functions in the proton [29]], 



u(x) = Z u u (x) + P un -/ d ® d (x) + K/7T+ ® P-K+d/u ® u (x) + -P U7T o /u g u (x) 

+ V u/K + g Pk+s/u <8> u (x) + ^V u/w o g [P n o u/u g u (x) + P^ d/d g d (x)] , 

d(x) = Z d d {x) + P d n+/u ® m (x) + V d/7r - g P^-u/d g do(aj) + -^Pd^/d ® dofa) 

1 

+ V^o g g d (a?) + -Vd/^ g [P ff o u/u g ito(a;) + P^d/d ® doOc)] , 

u(ar) = V^/ w - g Pv-u/d ® d (x) + ^V u/7T o g [P^/,, g w (^) + P ff o d/d g rf (a;)] , 

= V d/lT+ g P^+d/u g u (x) + -V s/7r0 g [P^o^ g u (x) + Pifidfd ® d (x)] , (16) 

where the constituent quark-distributions uq and do are normalized to two and one, respec- 
tively. Convolution integrals are defined as 



Pj a /i ®qi = I — Pj a /i {y)q%(-j 
j x y w j 

1 dyi f l dy 2 



V k / a g P aj /i g % -- 
In addition, V^./ a (x) follows the relationship 
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(18) 



We postulate that the bare-quark distributions are isospin-symmetric between the proton 
and the neutron, so we can obtain the quark distributions of the neutron by interchanging uq 
and do. Employing the quark distributions of the chiral quark model, we get the Gottfried 
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sum determined by the difference between the proton and the neutron structure functions, 

* da; 



Sg = f — [Fi(x) - F»(x)] 
Jo x 

l r 1 

= - dx [Au p (x) + Au p (x) - Au n (x) - Au n (x) + <P{x) + <P{x) - d n (x) - d n (x)] 
9 Jo 

= i + J dx / jj [u p (x) - u a {x)] + ~ [d*(x) - fr(x)] 



= l~l (Pr~) + I (P*+) + ^ «0> " <^°» • (19) 

We assume that the ISB is entirely from the mass difference between isospin multiplets. 
In this paper, we adopt {m u + m d )/2 = 330 MeV, m n ± = 139.6 MeV, m^o = 135 MeV, 
itik± = 493.7 MeV, and m^o = 497.6 MeV. We choose two sets of the mass difference 
between u and d quarks, namely 5m = 4 MeV and 5m = 8 MeV, respectively, in order to 
show the dependence on this important parameter. Based on Eq. ffTUl) and the experimental 
data of the Gottfried sum jr|, one can find that the appropriate value for is 1500 MeV. 
However, one cannot determine in the same method, because (Pk) in the Gottfried sum 
is canceled out. Usually, it is assumed that A^ = A^ = 1500 MeV 29|, |3J|. However, it is 
implied by the SU(3)j symmetry breaking that (Pk) should be smaller, and, accordingly, 
one should adopt a smaller Ak- In this paper, we adopt a wide range of A^ from 900 to 
1500 MeV. In addition, the parton distributions of mesons are the parametrization GRS98 



given by Gluck-Reya-Stratmann 32|, since the parametrization is more approximate to the 
actual value, 



V„(x) = 0.942aT a501 (l + 0.632^(1 - ^ a367 



V u/K+ {x) = V d/K o(x) = 0.541(1 - xf 17 V w (x). (20) 

We should point out that, in principle, it is possible that the parton distributions of different 
mesons in the same multiplet are different, and this can contribute to ISB simultaneously. 
However, in this paper, we simply neglect this possibility, and calculations in future can 
be improved if we have a better understanding of the quark structure of mesons. More- 
over, we have to specify constituent-quark distributions u and d , but there is no proper 
parametrization of them because they are not directly related to observable quantities in ex- 



periments. In this paper, we adopt the constituent-quark-model distributions [33] as inputs 
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for constituent-quark distributions. For the proton, we have 

2x Cl (l - x) Cl+C2+1 

UqIx) 



B[ci + 1,ci + c 2 + 2] 
x C2 (l - x) 2ci+l 



= B[ C2 + l,2c 1+ 2] - < 21 » 

where B[i,j] is the Euler beta function. Such distributions satisfy the number and the 
momentum sum rules 

i ,1 
u (x)dx = 2, / d (x)dx = 1, 



/ XMo(x)dx+ / xdo(^)dx = 1. (22) 
Jo Jo 



ci = 0.65 and c 2 = 0.35 are adopted in the calculation, following the original choice |33|, |34 |. 

We display the ISB of the valance- and sea-quark distributions in Figs, [fl El and |3j 
respectively. It is shown that in most regions, x5uy(x) > and x5u(x) > 0, and on the 
contrary that x5dy(x) < and x5d{x) < 0. Our predictions that x5u(x) > and x8d(x) < 
are consistent with the MRST parametrization 12 1| . and, moreover, the shapes of xSu(x) and 
x5d(x) are similar to the best phenomenological fitting results given by the MRST group. 
We should point out that our results are analogous to the results calculated in the framework 



of the meson cloudy model by Cao and Signal 



181 ]. and the shapes and magnitudes of x5u(x) 



and x5d{x) are similar to the results given in the framework of the radiatively generated 
ISB [191 ]. but with different signs. It can also be found that the difference between various 
choices of Ak is minor, but the different choices of 5m can have remarkable influence on 
the distributions. Especially, larger 5m can lead to larger ISB, and this is concordant with 
our principle that ISB results from the mass difference between isospin multiplets at both 
hadron and parton levels. From the figures, we can see that 5uy(x) reaches a maximum 
value at x ~ 0.5, and 5dy(x) has a minimum value at x ~ 0.4. It should also be noted that 
5qy(x) (q = u,d) must have at least one zero point due to the valance-quark- normalization 
conditions. We should also point out that at large x, 5uy/uy m —5dy/dy, and this implies 
that the magnitudes of the ISB for uy and dy are almost the same, but with opposite 
signs. Moreover, although both flavor asymmetry in the nucleon sea and the ISB between 
the proton and the neutron can lead to the violation of the Gottfried sum rule, the main 
contribution is from the flavor asymmetry in the framework of the chiral quark model. 
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FIG. 1: The ISB of the liy-quark distribution x8uy{x) versus x in the chiral quark model with 
different inputs. The red solid line is the result with 8m = 4 MeV and Ak = 1500 MeV as inputs. 
The blue dashed line is the result with 8m = 8 MeV and Ak = 1500 MeV as inputs. The green 
dotted line is the result with 8m = 4 MeV and Ak = 900 MeV as inputs. 



III. THE CONTRIBUTION FROM ISOSPIN SYMMETRY BREAKING TO THE 
NUTEV ANOMALY 

The measured sin 2 9w by the NuTeV Collaboration is closely related to the Paschos- 



Wolfenstein (PW) ratio |35| 



p- _ (°Nc) (°Nc) 1 -2n / 9 o\ 

R - (O-(^c) =2- Sm " w ' (23) 

where (o"nc) ^ s ^ ne neutral-current-inclusive cross section for a neutrino on an isoscalar 
target. If we take the ISB between the proton and the neutron into account, we obtain 
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FIG. 2: The ISB of the dy-quark distribution x5dy(x) versus x in the chiral quark model with 
different inputs. The red solid line is the result with 5m = 4 MeV and A# = 1500 MeV as inputs. 
The blue dashed line is the result with 5m = 8 MeV and Ak = 1500 MeV as inputs. The green 
dotted line is the result with 5m = 4 MeV and Ak = 900 MeV as inputs. 



where <5-Rp S w is the correction from the ISB to the PW ratio and takes the form 



r pISB 
°- rt PW 



sin 2 6 



w 



Jo x 


5uy(x) — 5dy(x) 


dx 




Uy(x) + dy(x) 


dx 



(25) 



with uy(x) and dy(x) standing for valance-quark distributions of the proton. We show 
the renormalization constant Z, the total momentum fraction of valance quarks Qy = 
Jq x [uy(x) + dy(x)] dx, and the correction of the ISB to the NuTeV anomaly A_Rp^, with 
different 5m and as inputs in Table HI It can be found that the ISB correction is of 
the order of magnitude of 10 -3 and is more significant with a larger 5m or Ak- Our result 
is consistent with the ranee -0.009 < AiZpg < +0.007, which is derived based on the 
parametrization given by the MRST group 21] . We should stress that the correction is 
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FIG. 3: The ISB of the sea-quark distributions x5q(x) versus x in the chiral quark model. The red 
solid line and the blue dashed line are the behaviors of x5u(x), with 5m = 4 MeV and 5m = 8 MeV, 
respectively. The green dotted line and the orange dash-dotted line are the behaviors of x5d(x), 
with 5m = 4 MeV and 5m = 8 MeV, respectively. 



remarkable, since the NuTeV anomaly can be totally removed if Ai?pw = —0.005, and, 
consequently, we should pay special attention to ISB in such problem. It is also worthwhile 
to point out that the correction is in an opposite direction to remove the NuTeV anomaly 
in the chiral quark model. Such a conclusion is the same as that given in the baryon-meson 



fluctuation model 



20(, but the value is one or 2 orders of magnitude larger. Our result of 



the ISB correction to the NuTeV anomaly differs from the results in Refs. 13, Il9 |. 



IV. SUMMARY 



In this paper, we discuss the ISB of the valance-quark and the sea-quark distributions 
between the proton and the neutron in the framework of the chiral quark model. We assume 



ff 



TABLE I: The renormalization constant, the total momentum fraction of valance quarks, and the 
correction of the ISB to the NuTeV anomaly in the chiral quark model. 



5m (MeV) A K (MeV) Z u Z d Q v Ai?p® 

4 900 0.7497 0.7463 0.8451 0.0008 

4 1200 0.7220 0.7185 0.8222 0.0008 

4 1500 0.6932 0.6896 0.7985 0.0009 

8 900 0.7515 0.7444 0.8455 0.0016 

8 1200 0.7239 0.7165 0.8227 0.0017 

8 1500 0.6953 0.6874 0.7990 0.0019 



that isospin symmetry breaking is the result of mass differences between isospin multiplets. 
Then, we analyze the effects of isospin symmetry breaking on the Gottfried sum rule and 
the NuTeV anomaly. We show that, although both flavor asymmetry in the nucleon sea 
and the ISB between the proton and the neutron can lead to the violation of the Gottfried 
sum rule, the main contribution is from the flavor asymmetry in the framework of the chiral 
quark model. It is remarkable that our results of ISB for both the valence-quark and sea- 
quark distributions are consistent with the MRST parametrization of the ISB of valance- 
and sea-quark distributions. Moreover, we find that the correction to the NuTeV anomaly 
is in an opposite direction, so the NuTeV anomaly cannot be removed by isospin symmetry 
breaking in the chiral quark model. However, its influence is remarkable and should be taken 
into careful consideration. Therefore, it is important to do more precision experiments and 
careful theoretical studies on isospin symmetry breaking. 
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